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ABSTRACT 

It is currently common to use spatially unresolved multi-filter broad-band 
photometry to determine the masses of individual stellar clusters (and hence 
the cluster mass function, CMF). I analyze the stochastic effects introduced by 
the sampling of the stellar initial mass function (SIMF) in the derivation of the 
individual masses and the CMF and I establish that such effects are the largest 
contributor to the observational uncertainties. An analytical solution, valid in 
the limit where uncertainties are small, is provided to establish the range of 
cluster masses over which the CMF slope can be obtained with a given accuracy. 
The validity of the analytical solution is extended to higher mass uncertainties 
using Monte Carlo simulations and the Gamma approximation. The value of the 
Poisson mass is calculated for a large range of ages and a variety of filters for 
solar-metallicity clusters measured with single-filter photometry. A method that 
uses the code CHORIZOS is presented to simultaneously derive masses, ages, 
and extinctions. The classical method of using unweighted UBV photometry 
to simultaneously establish ages and extinctions of stellar clusters is found to be 
unreliable for clusters older than k, 30 Ma, even for relatively large cluster masses. 
On the other hand, augmenting the filter set to include longer-wavelength filters 
and using weights for each filter increases the range of masses and ages that can be 
accurately measured with unresolved photometry. Nevertheless, a relatively large 
range of masses and ages is found to be dominated by SIMF sampling effects that 
render the observed masses useless, even when using UBVRIJHK photometry. 
A revision of some literature results affected by these effects is presented and 
possible solutions for future observations and analyses are suggested. 
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Introduction 



This paper is the third one of a series where we explore the effects of different biases on 
the determii iation of the stellar and cluster mass functions (SMFs and CMFs, respectively). 
In paper I (IMaiz Apellaniz fc Ubedal 120051 ) we analyzed the numerical biases induced by 
using bins of equal width when fitting power-laws to binned data (an effect that is more 
general than its application to the calculation of mass functions). Those biases can be 
eliminated in several ways, of which a simple one is b y grouping the data i n equal-number 
bins (as opposed to equal-width bins). In paper II (IMafz Apellanig |2008| ) I explored the 
effect of unresolved multiple systems, either physical or chance alignments, especially for the 
high-mass end of the stellar initial mass function (SIMF). In this paper I analyze the effect 
of random uncertainties in the mass determinations of individual stellar clusters and on the 
global properties of the obtained CMF as derived from spatially integrated (i.e. unresolved) 
photometry. I am currently working on the fourth paper of the series, which will explore the 
same issues as this one but referred to stellar instead of cluster masses. 

The measurement of the masses of unresolved stellar clusters in external galaxies has 
become popular in the last decade, especially thanks to the availability of HST imaging (see 
Zhang fc Fallll999l : lLarsenll2002l : Ide Griis et allboosl : lubeda et al.1l2007bl : IPowell et al II2OO8 
for examples). Accurately measuring stellar cluster masses is crucial to understand their 
evolution and, more specifically, their destruction rates and mechanisms. This is usually 
done by analyzing a large ensemble of clusters within a galaxy and deriving the present-day 
CMFs as a function of age. Current results regarding how and at what rate stellar clusters 
are destro yed are inconc lusive, with two alternative empirical models being proposed in the 
literature (iLamersI |2008| and references therein) . 



One outstanding issue with the calculation of stellar cl uster masses is SIMF samp hng. 
In a series of papers, Miguel Cervino and his collaborators (jCerviiio fc Luridianall2006l and 
references therein) have shown that for clusters with less than a certain number of stars 
the SIMF is not well sampled and, as a consequence, there are added uncertainties in the 
derivation of cluster properties from unresolved photometric or spectroscopic data. Thus, two 
clusters of the same mass, age, and metallicity can have quite different integrated properties 
because of the differences in their initial stellar population caused by the stochastic nature 
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of star formatior0. Undoubtedly, it is important to determine whether such effects are 
distorting the analysis of stellar cluster evolution by inducing biases in the observed CMFs. 

The papers by Cerviiio et al. deal mostly with the theoretical issues of SIMF sampling 
and provide a framework for its general study. The goal of this paper is more limited but at 
the same time more practical. It is more limited because I will concentrate on one observable, 
the cluster mass, and will deal with others (e.g. age) only as long as they influence the value 
of the observed mass. Furthermore, for the sake of simplicity I will restrict the analysis to 
solar-metallicity clusters. The practical character of this paper arises from its main aims, 
which are to provide [a] specific criteria to determine when cluster masses can be obtained, 
[b] methods for a more precise measurement of individual cluster masses, and [c] corrections 
to eliminate biases in both the individual masses and the CMF. In a sense, this paper 
can be understood as an unresolved-population equivalent to the resolved-stellar-population 
methods and tools used for extracting unbiased age s or star-formation histo r ies from color- 



magnitude diagrams developed in the last decade (Harris fc ZaritskyI 1200 ll : iDolphinI 12002 



J0rgensen fc Lindegrenll2005l : iNaylor fc Jeffries! l2006l ). 



The paper is organized in the following way. I start by briefiy describing the problem. 
Then, I select a simplified case (clusters with fixed single age, metallicity, and extinction) and 
derive an analytical solution for the effects of SIMF uncertainties on individual cluster masses 
and the CMF for the case where uncertainties are small and masses are derived from ^-band 
photometry (an appendix provides the mathematical formalism). Later, I use Monte Carlo 
simulations to analyze arbitrary mass uncertainties and develop an approximation using the 
Gamma distribution to obtain the behavior of the observed mass distribution for a given 
arbitrary cluster mass. The simulations are then extended to other photometric bands and 
known ages. Finally, the cases where the age or the age and the extinction are unknown and 
have to be derived from the same multi-band photometry as the mass are considered, and 
a method to obtain all of them simultaneously is presented and analyzed. I wrap up with a 
discussion on how these results apply to real data and objects and a list of conclusions. 



2. Description of the problem 

The standard technique to derive an observed cluster mass rrio with known age, ex- 
tinction, distance, and metallicity from spatially-integrated photometric data has two steps. 



^Note that other factors such as the dynamical evolution due to both internal and external causes can 
and indeed frequently do induce differences between otherwise initially similar clusters. However, those will 
not be considered in this paper. 
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First, the photometry is processed to account for distance and extinction to derive the lu- 
minosity in a given band (e.g. Johnson V), Ly- Second, evolutionary synthesis models are 
used to derive Qy-, the mass-to-luminosity ratio in that band for very large (real) cluster 
masses m, as a function of cluster age: 

Qy = hm — . (1) 
The observed cluster masses are then: 



rric 



QyLy. 



(2) 



What is the uncertainty a associated to the measurement of mo? There are three 
sources of uncertainty to consider: photometric, model, and SIMF sampling. For CCD mea- 
surements, the photometric uncertainty arises in most cases from the Poisson statistics of the 
detected radiation (but note that in some cases read noise and background subtraction can 
also play an effect). Under model uncertainties I include a number of effects that originate 
in our transformation of the photometry to Ly and the value of Qy, such as an inadequate 
characterization of evolutionary tracks and stellar atmospheres, incorrect calculation of ages 
and extinctions, and errors in the zero points of the photometric systems. The third uncer- 
tainty source, SIMF sampling, arises because the SIMF (understood as a probability density 
function) is only well sampled at very large masses; in other cases two clusters of the same 
mass can have significantly different populations of massive stars, which are the most scarce 
ones but can dominate the light output from young clusters. Therefore, SIMF sampling in 
a cluster dictates that two clusters of the same Ly, age, and metallicity do not necessarily 
have the same mass, an effect that becomes especially important for low-mass clusters. 

Which uncertainty source dominates? The ultimate answer will depend, of course, on 
the data and methods used to derive the observed masses. The behavior of model uncer- 
tainties is difficult to calculate because of the many parameters invo lved . However , given 



that photometric zero points are usually known to within 2% or better ( Maiz Apellanij 



20071 ). that current evolutionary synthesis models such as Starburst 99 Jleitherer et al.lll999h 



2006 



are quite capable of accurately reproducing the integrated spectra of massive clusters of dif- 
ferent ages, and that nowadays relativ ely sophisticated niethods can be used t o compare 



photometric data to model predictions ([Anders et al.ll2004l : iMaiz Apellanig|2004l ). it is safe 



to put a 5% cap on model uncertaintie 



;ions u 
intiecl. 



^If SIMF sampling is not relevant. As we shall see later on, some of the model uncertainties (e.g. 
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To characterize the other two types of uncertainties (photometric and SIMF samphng), 
I first consider a simphfied example that is relevant to current research: a family of single- 
age (10 Ma), solar metallicity, constant extincti on (Ay = 1.0) stellar clusters located at 
the distance of the Antennae galaxies (22 Mpc, ISchweizer et al.l |2008| ) observed with the 
Hubble Space Telescope ACS/WFC camera using the F550M filter with two exposures of 
1250 s each. Thephotometric uncertainties were calculated using the ACS Imaging Exposure 
Time Calculatocl. The SIMF sampling uncertainties were calc ulated using the evolu tionary 
synthesis add-on module available in version 3.1 of CHORIZOS (IMafz Apellanial2004j ) . which 
also produced the Monte Carlo simulations that are described later on. Results are shown in 
Table [T|, with both a low-mass and a high-mass cluster as examples, and a clear conclusion 
can be extracted from them: at least in this simplified cas^ SIMF sampling is the largest 
source of uncertainty in the determination of the masses of unresolved stellar 
clusters for masses lower than 10^ Mq and is clearly dominant in the more interesting (and 
currently actively debated) range 10^ — 10^ Mq. 



3. An analytical solution for /mo(^o) when a <^m 

As a first step towards the characterization of the possible biases caused by mass uncer- 
tainties in stellar cluster measurements, I develop an analytical solution that describes the 
behavior for the case a/m -C 1. I will consider that the real stellar cluster mass distribution, 
/m(^), is given by a power law with slope 7 (see Eqn. IA6I and Appendix A for the Bayesian 
formalism used in this paper) and that a has a power-law dependence on mass: 

(T(m) = a^(m/m,.)^, (3) 

where rrir is a reference mass, the uncertainty for that reference mass, and (3 a power-law 
exponent, with f3 < 1. We then have that the joint distribution of observed and real masses 
(see Eqn [M]) is: 



extinction) can be coupled with SIMF sampling issues and, thus, be greatly enhanced, 
• ^http : //etc . stsci . edu/webetc/acsI magingET C . j sp| 

*But likely also in many others, since for objects closer than the Antennae photometric uncertainties 
should be smaller. 
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A 



v27r(Tr(m/mr)^ 

A m^^mo^~^ 



exp 



mo — m 



m ' 



27r CTr 



2 V o"r(m/mi.)^^ 
exp [-aV(l + l/)-'''/2](l + Z/r^ 



(4) 
(5) 



where I have defined: 



y = {m - rrio) / rrio (6) 
a = m^^-^mj^ja,. (7) 

In order to guarantee that a/m <^ 1 I will consider only the range nio > with the 
additional assumption that ar/m^ -C 1. In that case, it is easy to see that a ^ 1 and that 
/mo,m(^o; '"^) is non-negligible only in those regions where \ay\ < 1, i.e. where \y\ ^ 1 and 
the observed mass is not too different from the real one. Under such circumstances we can 
do a Taylor expansion of the terms in (1 + to obtain: 



fMo,Mi^o,m) ^ = exp (-a y /2) x 



where I have kept all the terms up to Hence, the distribution of observed masses 

(see Eqn. I A3 II is given by: 

f r ^~ 7 A ^ (7 + 2^X7 + 2^-1) \ 

ImA'^o) ^ ( 1 + 1 . (9) 

I can now use Eqns. [S] and [3 to calculate first the mass correction for an individual 
cluster (c, see Eqn. lA4l) . and then the fractional number change (p, see Eqn. IA5I) and the slope 
change (A7, see Eqn. IA7p between the real and observed distributions for the approximation 
described above: 
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a" 



(10) 



V 



(7 + 2/3)(7 + 2/3 - 1) _ (7 + 2/3)(7 + 2/3-1) 
2a2 2 



(11) 



A. ^ _ (l-ffl7.2^)(7^2^-l) ^ _ , ,,)(^, , _ 1) (^^)^ (12) 

Equations [TUl [TT| and [T^ reveal that, to first order, the effect of uncertainties on c, p, 
and A7 depend on mo as mo"^*-^"^-* so that they tend to zero for large masses unless (3 = 1. 
The sign and specific magnitude of the effect of uncertainties on the three quantities depend 
on the specific values of 7 and /?. It is interesting to note that for 7 = —2/3, c, p, and A7 
are all zero, indicating that the observed and real distributions are identical and that no 
overall shift of the masses is observed to this degree of approximation. For 7 = 1 — 2/3, c > 
but p = A7 = 0, indicating that even though the two distributions are identical there is an 
overall shift from higher to lower masses for a given value i.e. real masses are larger than 
observed ones on averag^. For /3 = 1, Eqn. [12] indicates that the observed distribution is 
always a power law with the same slope as the real one and Eqns. [10] and [TT] give values of 
c and p independent of nio. 

When uncertainties originate from SIMF sampling (/3 = 0.5), Eqns. [T0| [TTl and [12] give: 



7 + 1 (7 + l)a. 



m^mo 



p 



7(7 + 1) 7(7 + 1)^^2 



2a2 



(13) 
(14) 



^^^_7(7 + l)^_7(7 + lK!_ (15) 

Therefore, in the low-uncertainty limit SIMF sampling generates mass corrections for 
individual clusters, fractional number changes, and CMF slope changes for large masses that 
are inversely proportional to the observed masses. For 7 < — 1, c and A7 are negative and 
p is positive. 



■^Of course, this is possible because I am not considering effects near the edges of fhii'm)- 
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4. Monte Carlo simulations for j\/(^o|^) 

The analytical solution developed in the previous section can help us characterize the 
effect of uncertainties on the observed mass distribution for stellar clusters but, considered 
by itself, it has two shortcomings. The first one is the scale of the effect e.g. if we rewrite 
Eqn.Ofor p = 0.5 as: 



a 



(16) 



we still need to calculate the value of 9 for the case of interest. For a Poisson process, Eqn. [in] 
is satisfied, so I will call 9 the Poisson mass. 9 can be thought of as characteristic cluster 
mass since for m = 9, a = m, and will be used extensively in this paper as a measurement of 
the stochasticity of clusters as a function of age and other parameters. The choice of notation 
will become clear when I introduce the Gamma approximation in the next subsection. The 
second shortcoming is that the analytical solution is expected to work only when the effect 
of uncertainties is relatively small and when we are using a power law without mass limits 
(which, of course, is unphysical because such distribution cannot be normalized). In order to 
evaluate its range of validity and to consider a more realistic case with possible large values 
of cr/m and mass limits, we need to go beyond an analytical approach and use Monte Carlo 
simulations of /mo|m(^o|^)- Such simulations can also provide us with the value of 9 and 
test whether Eqn. [16] is indeed correct. In this section and the following one I develop an 
example that, though a simplification of a real case, retains its most important effects. Later 
on I include some of the additional complications related to real data. 

Since I have already established the age, metallicity, extinction, and distance of our 
example in section [2l the next parameter that needs to be determined is 7. A number of 
indep endent results indicate that the CMF ha s a slope close to —2.0 in many circumstances 



Elson fc Falll Il985bl : IZhang fc Falll Il999l : lOev et al.l 120041 ) . so I fix 7 = -2.0 for the 



(e.g. 

simulations. For that value of 7 the analytical approximation gives c = — 1/a^, p = 1/a^ 



and A7 = — l/a^. I also fix m 



upper 



10^ Mq and let miower vary between 1 and 1000 M, 



0- 



The next issue to consider is the type of Monte Carlo simulation. Simula tions of stellar 



clust e rs typically use either a fix ed number of stars or a fixed total stellar mass (ICervino et al. 
2OO2I : IWeidner fc Kroupall2006l ). Here our goal is to generate /mo|m(^o|^), so I use Monte 



Carlo simulations of the second type. We also need to select the photometric band to 
calculate the observed mass (Eqn. [2]): our choice here is the most popular one, Johnson 
V. To generate the synthetic clusters I use t he evolutionary synth esis code available as an 
add-on to the latest version of CHORIZOS (IMaiz Apellanig 120041) . Each simulated 10 Ma 
cluster is a random realization of a Kroupa IMF (jKroupall200ll ) between 0.1 and 120 Mq 
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(note that it is possible that s ome stars may have aheady exploded as SNe) using a Padova 
isochrone (IMarigo et al.ll2008l ) for m < 7 Mq and a Geneva isochrone (iLejeune fc Schaerer 
2OOII ) for m > 7 Mq. The massive end of such an isochrone corresponds to a red supergiant 
with an initial mass of 18.29 Mq and My = —6.13. The observed mass associated with the 
lowest luminosity limit mo,LLL (i-e. th e maximum observed mass that can be produced by 
a "single-star cluster" of that age, see ICervino fc Luridianall2004l who use for such a 

quantity) for such a cluster is 2795 Mq, with the "single-star cluster" being in reality a 16.94 
Mq blue supergiant with My = —7.47. 

I generate Monte Carlo simulations of /a/o|m(^o|"^) for eight values of m between 3 ■ 10^ 
Mq and 1 ■ 10^ Mq. rig, the number of simulations, was adjusted for each of the eight values 
of m to guarantee a good sampling of /mo|m("^o|"^)- The output is shown in Tableland in 

Fig.m 

The first result is that ruo indeed behaves as a Poisson-like variable, with an average 
value fUo almost identical to m and a standard deviation a that satisfies Eqn. [16] to a high 
degree, as evidenced by the near-constant value of 6 derived for different masses. From 
Table [2] I obtain an average valu^ for 6 of 566 Mq, which will be used subsequently. 

The second result is that, for low values of m, fMo\M{^o\^) highly deviates from a Gaus- 
sian (s ee top plots of Fig. [1]) . Such a phenomenon has been described by lCervino fc Luridiana 
(120061 ) and it was easily predictable (but too often neglected in the literature): as m de- 
creases, a/m increases, with a = m for m = 6 and a > m for lower values of m. Under such 
circumstances a Gaussian cannot provide an accurate description because it produces a sig- 
nificant fraction of stars with negative masses, which is obviously unphysical. Furthermore, 
the complex shape of the luminosity function produced by the existence of fast post-MS 
evolutionary phases (e.g. yellow supergiants, see lCervifio fc Luridianall2006l ) induces the ex- 
istence of complex structures in /Afo|Af("^o|'^) for low values of m (e.g. the peak seen around 
1000 Mq in the upper left plot of Fig. [H which is real and not the result of using a small 
number for Ug). For higher real masses a Gaussian produces a reasonable approximation of 
fMo\M{fno\fn), as evidenced in the bottom plots of Fig. [H 

The non-Gaussianity of fMo\M{^o\^) introduces a problem. Eqns. lAll IA2[ and I A3 1 are 



no longer valid and since our goal is to calculate / 



Mo 



rric 



defined as: 



Mo 



fMo,Mi'^o,m)dm= / fMo\Mifno\m)fM{m) dm, 



(17) 



^This value, of course, is only applicable to the type of cluster and mass-measurement method described 
here. In other circumstances 9 should be different, as we shall see later on. 
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we need to be able to characterize fMo\M{^o\^) in a fine grid, i.e. to compute fMo\M{fno\m) 
for an arbitrary value of m within a reasonable amount of computing time. Monte Carlo 
simulations are computationally expensive (the time needed to generate the eight simulations 
in Table [2] was several hours) and generating hundreds or thousands of them may be simply 
prohibitive, so an alternative is needed. This is dealt with in the next section. 



5. The Gamma approximation and a solution for f'Moi'f^o) 



One approach to the generation of /Afo|Af('^o|'^) in a fine grid is that of lCervifio fc Luridiana 



(120061 ). Those authors suggest going beyond the Gaussian approximation by using an Edge- 
worth series that provides the correct values for the skewness and kurtosis of the distribution 
function. An alternative approach is to find a family of distributions that has the desirable 
properties of: (a) having the variance equal to the square of the mean; (b) being zero for 
negative values of the observed mass; (c) reasonably approximating the skewness, kurtosis, 
and asymptotic behavior for rrio and rrio ^ oo of the Monte Carlo simulations; and (d) 
behaving like a Gaussian in the limit m ^ oo. One choice is the Gamma distribution, a 
two-parameter distribution defined by the probability density function: 

-x/d 



fix;k,e) = X 



fc-i ^ 



where k and 9 are the scale and shape parameters (both positive), respectively, and T{k) is 
the normalization factor: 



r(A;) = / e-^e-'dt. (19) 

The Gamma distribution has mean kd and variance kO'^. Comparing this with the 
desired properties for /mo|m("^o|'^), I find that 6 has to be defined as in Eqn. [T6l and that 
k = m/6. Hence, the Gamma-like distribution function approximation for /mo1m(^o|'^) 
with the appropriate mean and variance can be written as: 

g-m,o/S 

/M„|M(mo|m) = ^o-/^-^ g^/,r(m/g) ' 



This expression is easily computable and should be understood as a function of mo, with 
m and 6 as parameters (i.e. it is a distribution of observed masses for given real masses and 
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Poisson mass 6). How good is this Gamma approximation? Figure [T] shows it for four dif- 
ferent values of m and Table [2] gives the values for k {= m/6, with 6 being the average value 
previously calculated), 71 (the skewness from the Monte Carlo simulations), 7i^r (the skew- 
ness from the Gamma approximation), 72 (the kurtosis from the Monte Carlo simulations), 
and 72,r (the kurtosis from the Gamma approximation). It can be seen that Eqn. [20] satis- 
fies the previously stated desirable properties, only underestimating the skewness by ~ 25% 
and the kurtosis by ~ 35%. Also, since the Gamma distribution is a smooth, single-peaked 
distribution it cannot reproduce the complex structures seen in the Monte Carlo simulations 
for low values of m but note that those peaks are not centered at the same values of rrio 
for different values of m (e.g. compare the region around 1000 Mq for the two top plots 
of Fig. [1]). Therefore, when one integrates over m in Eqn. [20] those small deviations from 
the gamma approximation should be mostly smoothed out. On the other hand, the asymp- 
totic behavior for small and large values of irio is excellent, and this is important in order 
to correctly characterize how many low-mass (e.g ~ 300 Mq) clusters are masquerading as 
intermediate-mass (e.g. 3000 Mq) clusters through sampling effects. 

Using the Gamma approximation for /mo|m(^o|^), I apply Eqn. [T7]to generate fMoif^o) 
for four values of miower between 1 and 1000 Mq (note that 1 Mq and 1000 Mq themselves are 
rather unrealistic extremes for miower, while 10 M0 and 100 Mq are more reasonable values). 
Results are presented in Fig. [2], which show the observed and real distributions (top) and 
the slope of the observed distributions along with the analytical solution. As expected, 
the analytical solution is a good approximation for nio ^ 10^ M©, since a/m = 0.24 for 
m = 10^ Mq, but fails for low- mass clusters. 

Figure [2] shows that fuoij^o) is relatively similar to fM{rn) for rrio > 10^ Mq (as expected 
from the analytical solution) but that there are significant differences for lower observed 
masses. First, in all cases except miower = 1000 Mq there is a "hump" (overdensity of 
/mo("^o) with respect to fMijn)) at ~ 1000 Mq which is generated by lower- mass clusters 
with relatively large apparent masses (due to the existence of one or several very bright 
stars). The height of the "hump" increases as miower decreases. Second, to the left of the 
hump fMoi^o) asymptotically tends to a power law with 7' > —1.0, with 7' decreasing with 
miower- Third, in order to connect the two asymptotic behaviors (7' > —1.0 for low observed 
masses and A7 = —O/irto for high observed masses) it is necessary for the observed slope 
7' to have a minimum with a value below —2.0 at a few thousand solar masses. As seen in 
Fig. [21 the depth of that minimum increases with decreasing miower- 



6. Filter, age, and extinction effects 



The case analyzed in the previous sections is quite simplified since we are using one single 
photometric band and we assume fixed known values for the age, metallicity, extinction, and 
distance to the clusters. In a more realistic case, one uses two or more filters and has to 
determine some or all of those four observables at the same time as the mass. It is easy 
to envision a case where all those quantities are variable but, in order not to excessively 
complicate the analysis, I will assume that the distance is known and the metallicity is 
fixed to solar, so the only possible variables are age and extinction. Also, only a uniform 
foreground screen far away from the cluster will be used to model the extinction. 



6.1. Other single filters and known ages 

In order to extend the analysis beyond the results for 10 Ma clusters using V derived 
in the previous section, I start by repeating the Monte Carlo simulations of /Afo|M(^o|^) 
for seven additional filters (Johnson U and B, Cousins R and /, and 2MASS J, H, and K) 
and eight additional ages, 1 Ma, 3.16 Ma, 31.6 Ma, 100 Ma, 316 Ma, 1 Ga, 3.16 Ga, and 10 
Ga (i.e. logarithmic values of 6.0, 6.5, 7.5, 8.0, 8.5, 9.0, 9.5, and 10.0 years, respectively), 
all of them with solar metallicity. The isochrones were derived from the same sources as 
the 10 Ma one and th e throughputs and zero points for the filters are those described in 



Maiz Apellanig (120071 ). The most luminous stars for the first three isochrones are main- 
sequence stars (MS, 1 Ma), luminous blue variables (LBV, 3.16 Ma)0, and red supergiants 
(RSG, 10 Ma), respectively, while for the last six (between 31.6 Ma and 10 Ga) they are 
asymptotic giant branch stars (AGB). Accordingly, I will refer to each of the nine ages as 
the MS, LBV, RSG, AGBl, AGB2, AGB3, AGB4, AGB5, and AGB6 phases, respectiveljfl. 



^Several points need t o be mentioned regarding the 3.16 Ma case. First, the isochrone pubhshed by 
Leieune fc Schaeren (j200l[ ) stops near the point where stars become LBVs and does not include more evolved 
(WR) stars. I address this point by interpolating the associated evolutionary tracks to extend the isochrone to 
the end of the Wolf-Rayet phase. Second, the evolution near the LBV stage is quite fast and, to some degree, 
unknown, with the star being capable of rapidly changing its surfac e temperature and sw itching between 
BSG, ASG, and WR spectral types even within a human lifetime (see Walborn et al. I I2OO8 I for an example). 
Third, it is not even clear whether all stars survive the LBV stage to become WRs or instead explode as 
SNe while they are LBVs (jSmith et al.l 120071 ). In any case, I verified that none of the above fundamentally 
affects the results in this section by running comparison simulations with the post-LBV phases excluded. 

^Note that the phases are named after the most luminous star possible for each age, not after the most 
luminous star that exists in a given cluster (i.e. a realization of the SIMF for that age). For example, a 
low-mass 10 Ma cluster can have only main-sequence stars because it does not have stars of the right mass 
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Figure [3] shows the isochrones for seven of the nine ages. 

For all filters and ages, I find that rrio behaves as a Poisson variable and that the Gamma 
distribution is a good approximation to /j\/o|a/('^o|^) over a large range of real masses. The 
values of 9 (Table [3]) show a large spread between 23 M© and 5717 M0, with two clear 
patterns: [1] As a function of age, 6 starts being small in the MS phase, grows very rapidly 
when stars start evolving and the cluster enters the LBV phase, reaches a maximum either 
in the LBV (optical bands) or RSG (NIR bands) phase, decreases until a mid-AGB phase 
and then increases again until the last age analyzed. [2] With the only exception of the MS 
phase, 6 increases as a a function of the effective wavelength, with the effect becoming larger 
with age. The trend in wavelength is reversed in the initial MS phase but the effect there is 
relative small {6k/0u = 0.6 at 1 Ma but 247 at 1 Ga). 

What is the reason for the behavior of 6* as a function of age and wavelength? One 
way to analyze the problem is to compute for a given isochrone and photometric band the 
product of the SIMF and the luminosity in that filter and to obtain the stellar mass at which 
the median of such a distribution takes place. That median luminosity point for a given 
filter tells us the location in the isochrone beyond which 50% of the light in that filter is 
produced. Then, to decide whether SIMF sampling effects are expected to be relevant in that 
band for a cluster with a given number of remaining stars we should compare that location 
to the point beyond which one expects only a small fraction of the stars to be present (i.e. 
the mass where the integral of the SIMF from there to the largest surviving mass is a small 
fraction of the integral for all the surviving masses). Such a comparison can be seen in Fig. [31 
which shows for the median luminosity points and the locations beyond where the fraction 
of surviving stars is 1/500 (0.2% of the stars) and 1/1000 (0.1% of the stars). 

The behavior of 6 is now easily understood. For the MS phase, the median luminosity 
points for all filters are clustered between the 1/500 and 1/1000 points, with the value 
increasing towards shortest wavelengths. Hence, I would expect relatively low values of 6, 
especially for the NIR bands. As we move to the LBV phase, two important changes take 
place: the most massive stars significantly increase their luminosity (with the rest of the 
stars remaining basically unchanged) and their surface temperatures switch from increasing 
to decreasing with mas^. This causes the median luminosity points to rise beyond the 
1/1000 point (producing a large increase in 6) and their order as a function of wavelength 
to reverse (inverting the tendency of 6 with wavelength). Nevertheless, in both the MS and 



(close to 18 M0) to have a RSG at that precise age. 

^Excluding the post-LBV phase, which contains a very low fraction of the stars and of generally lower 
luminosity. 
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the LBV phase the wavelength dependence is relatively small because all median luminosity 
points are clustered together. This changes in the RSG phase, where the large spread in 
temperatures for the brightest surviving stars causes the median luminosity points to fall in 
very different points of the isochrone. The U and B bands fall well before the 1/500 point 
(thus returning to the low values of the MS phase) while the R, I, and NIR bands remain at 
the right hand of the isochrone, well beyond the 1/1000 point. As we progress from 10 Ma to 
100 Ma, the median luminosities move left and down (i.e. towards less evolved stellar stages) 
with respect to equivalent points in the isochrone. This happens because the temperature 
difference between '"blue" and "red" stars decreases as a function of age while the luminosity 
difference between RG/AGB and MS stars is kept moderately low and causes all the values of 
6 to decrease as a function of age. For older clusters the trend in 6 is reversed because of the 
large luminosity difference between evolved and MS stars. The effect especially important 
in the NIR since in the optical it is partially alleviated by the displacement of the 1/500 and 
1/1000 points towards the beginning of the RG part of the isochrone. 

One practical problem with 6 is that Monte Carlo simulations are required to calculate 
it. Is there a way to avoid such simulations and to obtain 6 (or at least an estimate) directly 
from the isochrone and the SIMF? Two qua ntities that can be easily c alculated and that 
quan t ify the stochasticity of rrio are mo,LLL (jCervifio fc Luridianal 1200J) and ngg (IBuzzoni 



19891 : ICervino et al.ll2002l ). the latter defined as the effective fraction of stars contributing to 



the luminosity in a given band X: 



J^27 SIMF (m)Lx(m)c/m) 
''''' " /rirSIMF(m)L2,(m)rfm " ^^^^ 

From Eqn. [161 one expects 9 to be approximately m^^i,^/nes, where av is the average 
stellar mass for the SIMF (in our case, 0.71 Mq). I show in Fig. H] the correspondence 
between, 9, mo,LLL, and m^^ av/'^-efr for the 72 Monte Carlo simulations in this subsection. 
"^o.LLL is indeed correlated with 9 but with a quite significant spread. This can also be seen 
in Tables |3] and m where the NIR values of 9 for 3.16 Ma and 10 Ma are similar but those 
of mo,LLL differ by an order of magnitude. On the other hand, wi^, av/'T-eS is a much better 
predictor of the value of 9. For the 72 values in Fig. HI the ratio of m^,_av/''^eff to 9 has a mean 
of 1.04 and a standard deviation of 0.17. Therefore, if all that is needed is a rough estimate 
of 9, m^^s^r/nes seems to do a reasonable job. 

With the 9 values in hand and using the Gamma approximation for /mo|m('^o|'^); I can 
now apply Eqn. [T7]to generate /mo("^o) for clusters of different ages observed with different 
filters. In order to explore the full range of effects in /mo(^o) I show the two extreme cases 
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of 6 in Table [3] (23 and 5717 M©) in Fig. [5] for the same values of miower and mupper 
previously used. The overall structure seen in Fig. [5] is similar to that of the top panel of 
Fig. [2] but with some significant differences. First, the "hump" changes position, with its 
central position located close to 26. Second, for fixed values of miower and mupper the height 
of the "hump" increases with 9. Third, the combination of the dependences of the "hump" 
height with miower (previously noted) and with makes it disappear for low values of 6 not 
only for miower = 1000 Mq but also for miower = 100 M©. On the other hand, the "hump" 
is visible for all values of miower for = 5717 Mq. Finally, to the left of the "hump" the 
value of 7' decreases as 6 increases. In all cases I find that the analytical solution provides 
the correct asymptotic solution for large values. Therefore, if one wants to fix a value for 
A7 (<^ 1) as a tolerance limit for the effects of SIMF, it is possible to use Eqns. l3]and[T6]to 
obtain that the observed masses have to satisfy: 

For example, if we want to tolerate slope changes up to 0.1 and 7 = —2, then mo > 10 6. 



6.2. Unknown ages 



So far, I have considered that the cluster ages were known prior to the measurement of 
their masses. In most cases, this is unrealistic since one typically uses the photometry to 
simultaneously determine mass and age. This is done by obtaining multi-band photometry, 
using the colors to first determine the age (and maybe also the extinction, which will be 
assumed to be non-existent in this subsection), and then using one of the ma gnitudes to 
calculate the cluster mass. A classical example is that of lElson fc Fall! (Il985al ). who used 
spatially-integrated UBV photometry to assign ages to a collection of LMC clusters by 
plotting their location in a. U — B versus B — V diagram. 

Working with color-color diagrams (e.g. U — B versus B — V) to calculate ages has the 
limitation that only two photometric quantities can be used simultaneously. In order to take 
full advantage of the p hotometric informatio n, one can use a multi-band mini mization code 
such as CHORIZOS jMaiz ApellanizI booi ) or AnalySED jAnders et al.ll2004[ ). Such codes 
allow one to explore a multi-dimensional parameter space while using an arbitrarily large 
number of photometric bands. Increasing the number of filters has the advantage of using 
additional information that can be especially useful when the observed photometry is similar 
to the model spectral energy distributions (SEDs) but not identical, such as when SIMF 
sampling effects are present. Thus, the additional filters can in some cases (at least partially) 
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"correct" the effect of one magnitude that is offset by such effects and provide a solution 
which is closer to the real one. Another advantage of a multi-band minimization code are 
the possibilities of assigning different weights to each filter and of estimating uncertainties 
in the output parameters. 

In this paper I will use CHORIZOS, a code i n which several change s have been included 
(its latest version is 3.1) since it was described in lMaiz ApellanizI (120041 ). CHORIZOS is now 
a full Bayesian code that allows the use of priors and that can work with either magnitudes, 
colors, or spectrophotometry. It also now has a unified atmosphere grid that combines recent 
SEDs from diff erent authors. More specifica lly, for O and B stars it uses the TLUSTY model 



anz &: Hubenvl (l2003l . 120071 ) and for late- type supergiants the MARCS model 



PlezI (|2008|). Version 3.1 of CHORIZOS has the added flexibihty of dealing 



atmospheres of 
atmospheres of 

with almost any combination of up to 5 parameters and has an evolutionary synthesis add-on 
module. 

The Monte Carlo simulations for the 10*^ M© clusters with ages between 1 Ma and 10 Ga 
previously described were used as input for CHORIZOS. Three types of executions depending 
on the photometry and weights used were performed: [a] UBVRIJHK photometry with 
weights for each filter deriv ed from the age-dependent single-band values of 9 in Table [3] (see 
Cervino fc Luridianal 120071 ) : [b] UBVRIJHK photometry with uniform weights for each 
filter; [c] UBV photometry with weights derived as in [a]. CHORIZOS was executed for 
the 10 000 simulated clusters generated for each real age and execution type by comparing 
the observed colors with the expected photometry of well-populated (i.e. of infinite mass) 
clusters with ages between 1 Ma and 10 Ga and evaluating the likelihood in each case. Such 
an experiment is a rather realistic simulation of what happens with real data. Note that our 
solution space is one-dimensional (age) while our data has 2 or 7 free inputs (colors) and that 
SIMF sampling effects will force the observed photometry to be incompatible in a strict sense 
with the available solutions. For each simulated cluster, the observed age was taken to be the 
mode of the likelihood distribution and the combination of the V magnitude and the mode 
of the observed age was used to derive its observed mass (taking into account the expected 
evolution of V with age for clusters without SIMF sampling effects). The distribution of 
observed ages (modes) for each real age and execution type is shown in Fig. O Table O gives 
the median of the observed age and mass distributions, as well as their inferior and superior 
uncertainties derived from their respective distributions. The last column in Table [5] gives 
the uncertainty in mass derived from the value of 6 for V in the previous subsection (i.e. 
assuming that age and extinction are known). 

I first compare executions [a] and [b]. For MS and LBV clusters, results are nearly 
undistinguishable. However, for later ages, the weighted executions are significantly better. 
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Unweighted executions have in general complex distributions in age, with multiple peaks 
present in most cases. Some of the secondary peaks contain more than 20% of the clusters 
and are located at observed ages very different from the real ones (e.g. log(age) of 6.8 
for 1 Ga clusters). The weighted executions with UBVRIJHK photometry, on the other 
hand have, for the most part, three desirable qualities: they are [a] nearly- Gaussian with 
[b] relatively narrow widths, and [c] centered close to the real age. The three (relatively 



(non-Gaussian distribution with moderately large width), and 100 Ma (secondary minor 
peak around log(age) of 6.6). Those three ages will be analyzed in the discussion but let 
us note here that the unweighted executions show the same or worse problems at those 
ages. Therefore, I conclude that using photometric weights based on the degree of 
stochasticity of each filter can significantly reduce the errors in the estimated 
ages of clusters. 

I now compare executions [a] and [c]. In general, differences between them are smaller 
than between [a] and [b] and can be broken down into three age regimes: 

1. For MS and LBV clusters, differences are minor, as it was also the case between [a] and 
[b] . The only significant difference is the presence of a few LBV clusters with observed 
log(ages) between 7.0 and 8.0 for the UBV photometry execution (barely visible at 
the bottom of the lower plot in Fig. [6]). 

2. For RSG and AGBl clusters, UBVRIJHK photometry has the advantage with respect 
to the determination of ages. For 10 Ma clusters, the UBV execution shows a strong 
secondary peak at log(age) of 6.6 that is much weaker in the UBVRIJHK execution. 
For AGBl clusters, both executions are relatively poor but the UBVRIJHK case has 
a somewhat wider distribution. 

3. For ages between 100 Ma and 10 Ga, the UBV executions provide better results with 
narrower widths in the distribution. 

Regarding the measured masses, two results can be derived from Table [51 First, the 
median masses for executions [a] and [c] are in most cases similar to each other and a few 
percent lower than the real masses. The most significant differences are for MS clusters (in 
both executions the median masses deviate from the real ones by ~13% and the distribution 



^°One point to be considered for 1 Ma and 10 Ga clusters is that they are located at the edges of the 
possible age ranges, thus possible offsetting the observed distribution by creating "one-sided Gaussians". 
However, as discussed later on, this is likely not the source of the offset for 1 Ma clusters. 



minor) exceptions are 1 Ma clusters 
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has a long tail toward even lower masses), RSG clusters (the median mass from [a] is much 
better than the one from [c]), and AGBl clusters (the median mass from [c] is much better 
than the one from [a]}. Second, the introduction of an unknown age does not affect much the 
dispersion in masse^lj for 3.16 Ma and 10 Ma clusters (compared to the dispersion using V 
photometry for cluster of known age) but for the rest of the ages it increases the dispersion 
by factors between 2 and 3. This single-mass analysis suggests that the value of 9 is left 
unchanged for LBV and RSG clusters but can increase up to an order of magnitude for 
clusters of other ages (for a fixed mass, 6 is proportional to o"^). In order to verify such 
an assumption, we need to analyze the equivalent results to those in Table [5] for a different 
mass. For that purpose, I have repeated the CHORIZOS analysis using the Monte Carlo 
simulations for 10^ M0 clusters and verified that increasing the cluster mass by a factor of 
ten reduces the dispersion in mass by factors between 3 and 4 (i.e. similar to the expected 
VTO), so it appears that, at least in this mass rang the uncertainties in cluster masses 
with unknown ages approximately satisfy Eqn [161 

One preliminary conclusion from this comparison between executions [a] and [c] is that 
even though adding more filters to the mix has potential benefits, this is not always the case. 
I will revisit this issue when extinction is included. 



6.3. Unknown ages and extinctions 

As a final step towards simulating a realistic measurement of masses for an ensemble 
of stellar clusters, I consider the case in which mass, age, and extinction are simultaneously 
determined from the photometry. Once again, I use the Monte Carlo simulations for 10^ Mq 
and 10^ Mq clusters. Having previously determined that photometric weights based on 
the degree of stochasticity significantly reduce the errors, I consider only the two types 
of CHORIZOS execution in the previous subsection with non-uniform weights (based on 
UBVRIJHK and UBV photometry, respectively). I extinguished the phot ometry for all 
the Monte Carlo simulations using a color exces ^£^(4405 - 5495) = 1.0 and a lCardelh et all 



^^In TableOl give inferior and superior uncertainties derived from percentiles instead of a single value for 
the dispersion (i.e. the standard deviation of the distribution) because some cases show highly asymmetric 
distributions. Nevertheless, in most cases the standard deviation is approximately the average of the two 
uncertainties given. 

i2por very low cluster masses, this is likely to break down because one would expect all distributions in 
age to become very wide and with multiple peaks. 

^•^£^(4405 — 5495) is the monochromatic equivalent to E{B — V) and is a direct linear measurement of the 
amount of dust for a given extinction law, as opposed to E{B — V)^ which is also a function of the input 
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( 119891 ) extinction law with -R5495 = 3.1. I executed CHORIZOS in an analogous way to the 
way it was done in the previous section but this time with two free parameters, log(age) 
between 6.0 and 10.0 and £^(4405 — 5495) between 0.0 and 2.0. The mode of the output 2-D 
likelihood table as used to select the age and color excess of each cluster. From those two 
values and the V magnitude the mass was finally calculated. The results for the ages are 
shown in graphical form in Figs. [7] and [8] and for the ages, color excesses, and masses for the 
10^ M0 case in tabular form in Table [61 

As opposed to the case where extinction is known, I find that, for most ages, the 
UBVRIJHK results are significantly better than those using UBV photometry. UBV 
photometry does a particularly bad job of determining the ages for solar metallicity clusters 
with ages of 31.6 Ma, 316 Ma, 1 Ga, and 3.16 Ga. In those cases multiple peaks are present 
in the age distribution and the observed masses can differ from real ones by large factors. 
For 100 Ma and 10 Ga clusters UBV photometry does a better job (even slightly better 
than UBVRIJHK) but for the youngest clusters (1 Ma to 10 Ma) it is significantly worse 
than UBVRIJHK, since the lower plots of Figs. [7] and [8] are the only cases where MS/LBV 
clusters are misidentified as RSG clusters and vice-versa. Therefore, I conclude that the 
addition of RIJHK photometry can significantly reduce errors in mass and age if 
extinction is unknown. 

The results for clusters between 1 Ma and 100 Ma using UBVRIJHK photometry in 
Table E] are very good: the uncertainties in mass are similar to the reference ones (where both 
age and extinction are known) and only the youngest clusters are slightly biased towards older 
ages and lighter masses (an effect that we already encountered in the previous subsection). 
For 100 Ma and 10 Ga clusters with 10"^ M©, UBVRIJHK still provides relatively unbiased 
ages and masses with mass uncertainties only a factor of two larger than the reference ones. 
The results of lower quality appear for clusters of 31.6 Ma, 316 Ma, and especially, 1 Ga 
and 3.16 Ga. For 31.6 Ma clusters of 10'^ Mq, the age distribution is three times wider than 
the reference one while for 316 Ma clusters of the same mass a secondary peak is present 
around 5 Ga. For 1 Ga and 3.16 Ga clusters of 10^ M0, the age distribution is very wide 
and biased; therefore, the age and mass cannot be accurately determined. The situation 
is significantly improved for 10^ M© clusters (Fig. [H]), where mostly unbiased and narrow 
distributions are observed for all ages using UBVRIJHK photometry (but not so for the 
UBV case). Therefore, if extinction is not independently known, t/5V7?/Ji7i^ photometry 
can be used to accurately measure ages and masses of unresolved stellar clusters only for 
some of the interesting regions of the age- mass 2-D space. 



SED and has a non- linear behavior. Similarly, i?5495 is the monochromatic equivalent to Ry- 
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7. Discussion 



In this paper we have seen that SIMF samphng has large effects in the measurement 
of stellar cluster masses and the associated CMFs which, if ignored, can lead to substantial 
biases on the derived results. Here I explore under which circumstances it is possible to 
extract valid and unbiased information about cluster masses. To that purpose, this discussion 
is divided into four parts: [a] a description of the cluster mass and age ranges over which 
masses can be accurately measured, [b] an exploration of some effects that have not been 
taken into account in previous sections, [c] a review of the validity of some of the literature 
results, and [d] guidelines for future work that may extend the measurement ranges in mass 
and age. 



7.1. What cluster masses and CMFs can be measured with broad-band 

photometry? 

The two most important results of the previous section are the need to give weights to 
different filters depending on their degree of stochasticity and the advantage of including RI 
and NIR filters besides the traditional UBV set when extinction is unknown. Therefore, for 
this subsection I will assume that UBVRIJHK photometry is used in combination with 
CHORIZOS or a similar code. 

The most encouraging results are for the youngest clusters, those between 1 Ma and 10 
Ma. The values of 9 lie between 500 Mq and 1000 Mq and are relatively impervious to the 
inclusion of extinction. Furthermore, even though some mixing is present between 1 Ma and 
3.16 Ma clusters for 10^ M© clusters, UBVRIJHK photometry clearly isolates them from 
older clusters. The mixing between the two youngest ages takes place (for 10^ M©) with a 
significant fraction of 1 Ma clusters and a much smaller fraction of 3.16 Ma clusters having 
observed ages around 2 Ma. The reason for the displacement of some the youngest clusters 
towards slightly older observed ages lies in the relative dearth (due to SIMF sampling effects) 
of very massive stars, which makes those clusters have U — B colors similar to those of 2 Ma 
cluster^ (the 2 Ma isochrone is not shown in Fig. [3] but is very similar for most masses to 
the 1 Ma isochrone, with the only difference taking place near the top of the HRD, where 
it is slightly displaced towards lower temperatures). Since a well-populated 2 Ma cluster 



^"^Note that there is one technique that is widely used by itself or in combination with broad-band pho- 
tometry to distinguish ages of young clusters and that is based on the use of nebular lines to estimate th e 
ionizing flux of the cluster and, from there, its age. See Maiz Apellaniz ( 2000l ): MacKentv et al. ( 2000 ): 
Fall et all (|2005h for examples. 
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is brighter than a 1 Ma cluster of the same mass in V (fundamentally due to the strong 
dependence of the bolometric correction on temperature), the observed masses for 1 Ma 
clusters are slightly biased towards lower values in Tables O and O 

Clusters in the LBV phase show an interesting characteristic: they have high values of 
6 for all filters in Table [3] but a relatively low value of the same quantity when UBVRIJHK 
photometry is used, even when extinction is allowed to vary. The explanation is threefold: 
[a] the brightest stars in such clusters (mostly B supergiants) have colors which are similar to 
the integrated cluster colors, so adding or subtracting a single bright star changes all of the 
magnitudes in a similar way; [b] the non-extinguished colors of 3.16 Ma clusters are rather 
unique (only clusters of 1 Ma have relatively similar colors), so it is hard to mistake their 
ages; and [c] extinction changes their colors in an also unique way. With those properties, 
the combination of properly-weighted multiple filters can yield a combined value of 9 lower 
than the one derived from a single filter such as V . 

RSG clusters share some of the properties of LBV clusters. In particular, their colors 
are unique because they are "blue in the blue" and "red in the red" i.e. their SEDs are 
dominated by early-B type stars in the optical and by late-type supergiants in the NIR (see 
Maiz Apellaniz et al.ll2004al for a practical application of this property). This makes 10^ M© 



clusters of ages close to 10 Ma relatively easy to distinguish from clusters of other ages, even 
after SIMF sampling effects and extinction are included. In summary, I conclude that the 
effects of SIMF sampling on the observed ages and masses of young clusters (1 Ma 
to 10 Ma) determined from broad-band photometry can be relatively minor. We 
can use Eqn. |22] to quantify this and specify that SIMF sampling effects introduce a change 
in the CMF slope of at most 0.1 for clusters in the MS. LBV, and RSG phases for masses of 
10^ Mq and higher if UBVRIJHK photometry is used with the appropriate weights. 

The situation changes as we leave the RSG phase. Age uncertainties for AGBl clusters 
are several times larger than for 10 Ma clusters of the same age and in most simulations the 
age distribution has multiple peaks. The decrease in precision is also noticed in the observed 
mass, though to a lower degree. Uncertainties decrease as we reach an age of 100 Ma, with 
values for the age uncertainties slightly worse than for young clusters and values for the 
mass uncertainties which are similar or even better. For clusters older than 100 Ma, the 
critical factor is whether extinction is independently known or not. In the former case, age 
uncertainties can be kept under control while in the latter they are very large already for 10^ 
Mq clusters. In summary, the effects of SIMF sampling on the observed ages and 
masses of intermediate-age and old clusters strongly depends on the specific age, 
filters used, and whether extinction is independently known or not, so a specific 
analysis is needed for each observational setup and circumstances. In some cases it will be 
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possible to keep the uncertainties under control but in others broad-band photometry will 
be almost useless to determine ages and masses. More specifically, post-RSG clusters 
(AGBl phase, log(age) ~ 7.5 Ma) are especially problematic. 

Even though it is possible to measure ages and masses for clusters in the MS, LBV, and 
RSG phase, one still has to be careful with contamination from older clusters. That can 
happen either because older ages have broad and smooth observed age distributions that 
produce an overlap in age or because of the existence of well-defined secondary peaks at the 
wrong age (e.g. 31.6 and 100 Ma clusters at log(age) of 6.6 and 31.6 Ma clusters at log(age) 
of 7.1 in t he top panel of Fig. [3 These secondary peaks can be observed when analyzing real 
data (see iBik et al.ll2003l for M51), though their exact location depends on the used models, 
because the beginning and end of the RSG phase are a strong function of the metallicity, 
rotation, and input physics used. The peaks take place at ages where colors change rapidly 
in time, thus creating large unique regions in color space where SIMF sampling effects can 
move clusters from other ages. The concentration at some ages can also create gaps in others, 
such as the one observed around log(age) of 7.1 in iFall et al.l (120051 ). Those authors indicate 
that the gap is caused by the introduction of observational errors in the fitting procedure 
but fail to notice the specific nature of those errors (SIMF sampling). 



7.2. Other possible problems with the CMF 

As we have seen, SIMF sampling affects not only the observed masses of individual 
clusters but also the CMF. I briefly describe here four effects that, when coupled with SIMF 
sampling, can introduce additional biases in the CMF. All of them have to be considered 
when deriving cluster mass functions from real data. 

The flrst effect is the magnitude limit. SIMF sampling complicates the derivation of 
completeness corrections because there is no longer a one-to-one correspondence between 
luminosity and mass for a given age. In general, completeness corrections are expected to 
be smoother functions due to this issue. 

The sec ond effect is the distinc tion between point sources and extended objects. Some 



authors (e.g. IWhitmore et al.lll999[ ) use radial concentration indices to separate real clusters 
from stars. However, such algorithms may identify a cluster dominated by a single star (e.g. 
a ~ 10^ Mq, 3.16 Ma cluster with the "right" sampling of the SIMF) as a point source, hence 
biasing the derived CMF. Later on I describe one possible way to address this issue. 



The third effect is the possible existence of multiple stellar generations within a cluster 
or association, a characteristic that has been observed from very young objects (e.g. 30 
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Dor, IWalborn et al.l |2002| ) to old ones (IMackey et al.l |2008| ) and that invalidates the use of 
single-burst models (i.e. isochrones) or at least degrades the quality of its results. Such 
multiple generations a re more likely to happ en in associations or massive compact clusters 
surrounded by a halo (IMaiz Apellanizll2001bl ). Those objects should be more easily resolved 
than real, compact clusters, which is one reason for the need of high spatial resolution to 
increase the accuracy of CMF studies. Some studies do not take into account that not all 
apparent clusters in distant galaxies are real clust ers but are instead associations which may 
contain a significant fraction of the massive stars (jGarmanylll994l : IChu fc Gruendlll2008l ) and 
that should dissolve rather rapidly, h ence providing a simple explanation for the phenomenon 
of infant mortality (IFall et al.ll2005l ). 



The fourth effect is the evolution of the SMF in a cluster with time due todynamical 
effects, whih invalidates equating it with the SIMP, even after taking into account stellar 
mass loss and the existence of stellar remnants. Nowadays it is recognized that there is a 
continuum in the total (kinetic + potential) energy of stellar groups at birth that spans 
the z ero energy point due to the dominating role of turbulence in the star formation pro- 
cess (iMac Low fc KlessenI I2004J : IClark et al.ll2005l ). Furthermore, the cluster m ay become 
unbound soon after formation due to gas expulsion (IGoodwin fc BastianI 120061) or later on 
by intra-cluster stellar encounters or external tidal effects ( Fall fc ReesI 19771 ). The loss of 
stars prior to comp lete dispersal can naodify the observed photometry significantly, espe- 
cially for old ages ( IKruijssen fc Lamerd l2008l ). The third (multiple stellar populations in 
clusters/associations) and the fourth (dynamical evolution) effects are tied together, as evi- 
denced by the difficulty in analyzing the bound character of some stellar groups. 

The fifth effect tied up with SIMF sampling is the validity of the stellar initial mass 
function itself and of the isochrones used for the comparison between the observed data 
and the model predictions. Since most of the luminosity is produced by the most massive 
surviving stars in a cluster, the major effect o f the discrepancie s between the real SIMF and 
the assumed one (in this paper, Kroupa, see lElmegreenl |2009| for a recent review on SIMF 
variations) should be a normalizing factor for the total mass, wh ich can be analyzed with 
veloc ity dispersion studies through the use of the virial theorem (IKouwenhoven &: de Grijs 
20081 and references therein). The validity of the isochrones is a more complex issue. For the 
youngest clusters, the two main possible sources of biases are [a] the interplay between rota- 



tion, metallicity, and mass loss in the post-MS evolution as a function of mass (IMeynet et al. 



20081 ): and [b] the rejuvenation of the SED induced by mass transfer in close binary systems 
that allows for the e xistence of O and/or WR stars beyond the 6 Ma age expected from 
single-star evolution ( van Bever et al. 1999 : Mas- Hesse &: Cervifiolll999l ). However, any of 
those two issues is unlikely to change the measured ages (either with broad-band photometry 
or with the equivalent width of Balmer lines) by more than a few Ma and should have a 
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relatively small effect in the measured masses. Older clusters pose most significant problems. 
First, the short-lived post-AGB phase (not included in the isochrones here) is sufficient to 
produce large increases in 6 for the short- wavelength filters because just an individual pAGB 
star can outshine the rest of the cluster. Second, at low metallicities horizontal-branch stars 
have a similar or larger effe ct at short wavelen gths (and if one goes into the UV, even white 
dwarfs come into play, see iKnigge et al.l l2008l for an example). Third, blue stragglers (the 
old-age equivalent of the rejuvenated O and WR stars previously mentioned) produced by 
collisions and binary mergers also alter significantly the short wavelength filters and their 
stochastic behavior. Blue stragglers have the additional difficulty of strongly depending on 
dynamical effects, thus complicating the construction of isochrones to model them. In sum- 
mary, the low values of 9 for U and B in Table [3] for old clusters are not realistic in most 
circumstances. 



7.3. Testing the validity of literature results 

Several works have been published in the last decade attempting to measure the CMF 
and its evolution with time. The results in this paper can be used for a qualitative exami- 
nation of their validity. 1 discus^ three of them here, all of them based on HST/WFPC2 
photometry: 



Chandar et aP Jl999al ffi used [/-BV^-equivalent + FUV filters to study the stellar 



clusters in M33 {d = 0.84 Mpc). 



Bik et al.l (120031 ) (see also iBoutloukos fc Lamersll2003l ) used ?7-BV^-RJ-equivalent filters 



to study the stellar clusters in M51 {d = 8.4 Mpc). 

Fall et al.l J2005h (see also Iwhitmore et~aD Il999l : Izhang fc Fall Il999h used UBVI- 



equivalent -|- Ho; filters to study the stellar clusters in the Antennae {d = 22 Mpc). 



A first criticism to all of the three analyses is that none of them uses weights based on 
the expected behavior of 6. This issue affects especially the M51 and Antennae papers, since 
they use a combination of filters that range from [/-like (F336W) to J-like (F814W), which 
differ greatly in that aspect. 



^■^Only strict SIMF sampling effects plus extinction are dealt with in this subsection; other issues such as 
distance, model atmospheres, photometric system conversions, metallicity effects, or analysis code differences 
are ignored. 
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The observed log(age) for the clusters analyzed by Chandar et al. in M33 range from 
less than 6.6 to 10.3, similar to the ages used in this paper. Also, most of their analysis 
is based on UBV photometry due to the low detection rate of objects in their sample in 
the FUV (due to the low sensitivity in that filter and the intrinsic weakness of the FUV 
luminosity for o ld clust ers) . Therefore, our third execution type is the relevant comparison. 



Chandar et al.l (jl999bl ) use external extinction determinations based on nearby stars and 



the overall extinction (both foreground and internal to M33) is low, so reddening corrections 
should be pretty accurate. They find that the majority of their sample has observed masses 
in the range 4 ■ 10^—10^ M0 and log(age) lower than 8.5. Hence, our results in Table [5] 
and Fig. [6] are lower limits on the uncertainties induced by SIMF sampling for the M33 
clus ters. It can be s een tha t SIMF- related uncertainties are relatively large and should bias 
the IChandar et al.l (Il999bl ) results. In particular, one would expect that clusters with real 
log(age) close to 7.5 would suffer especially from SIMF sampling effects and have different 
ob served ages from their r eal ones, thus biasing the observed age distribution. Indeed, Fig. 10 



m 



Chandar et al.l (jl999bl ) shows a minimum close to that age and an overall maximum in 



log(age) between 7.8 and 8.4, which points in that direction. 



Most of the M51 clusters analyzed by iBik et al.l (120031 ) have observed ages below 1 Ga, 
with less than 1% of them being apparently older. Their sample is significantly larger than 
the M33 one and some of their clusters are more massive but still more than 80% have 
masses below 10^ Mq. They use UBVRI photometry without stochasticity-based weights 
and simultaneously fit masses, ages, and extinctions. For the latter they obtain values of 
E[B — V) between 0.0 and 1.0. None of execution types in this paper fits exactly into 
this description since the ones where extinction is left as a free parameter use either UBV 
or UBVRIJHK photometry with stochasticity-based weights and the one where uniform 
weights are used is for UBVRIJHK photometry with known extinction. Nevertheless, 
we can use the following argument: [a] Since the use of uniform weights enhances the 
uncertainties associated with SIMF sampling, our unknown extinction executions (which 
use stoch a sticity -based weights) should provide an optimistic limit for the validity of the 



Bik et al.l (120031 ) results, [b] I would expect UBVRI photometry to provide an intermediate 



case between UBV and UBVRIJHK photometry, so analyzing both executions should give 
us an insight on the M51 results. With those ideas in mind, I first find that clusters with real 
log(age) near 7.5 should once again (as for M33) be scattered in observed log(age) between 
6.6 and 8.2. Fig. 13 indeed shows a local minimum between log(age) of 7.5 and 8.0. Also, as 
noted by the authors, local maxima are found around observed log(age) of 6.70, 7.45, and 
maybe 7.20. The specific locations of the maxima depend on the code used to generate the 
synthetic photometry and, therefore, they should not be overanalyzed. Nevertheless, their 
existence points towards the existence of SIMF sampling effects in the data. Furthermore, 
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this effect cannot be ignored as being of little importance, since, as we have seen, some of 
the clusters observed at those ages can have real ages which differ by ~1.0 in log(age). The 
situation for clusters older than 100 Ma is not good either. Even though their average mass 
is larger. Figs. [7] and [H] show us that when extinction is unknown, it is difficult to measure 
masses and ages of old clusters from integrated broad-band photometry even if their masses 
are as large as 10^ M0. 



The masses of the Antennae clusters in iFall et al.l (120051 ) are significantly higher than the 
ones in the previous studies, since they analyze clusters of more than 3 ■ 10^ Mq. They also 
incorporate Ha observations to their UBVI-like photometry which helps in the separation of 
clusters younger than log(age) lower than 6.6 from older ones. Those characteristics are likely 
to reduce the effects of SIMP sampling with respect to the M33 and M51 data. Nevertheless, 
it is clear from their Fig. 1 (as the authors acknowledge) that the observed age distribution 
has a nonrandom error component, with a strong concentration of clusters below log(age) of 
7.0 and a dearth between 7.0 and 7.4. As previously d iscussed, such a distribution is a clear 
sign of SIMP sampling effects and a look at our Fig. [8] (IFall et al.ll2005l fit the extinction from 
the same data used to calculate ages and masses) show that even for 10^ M© clusters, broad 
band photometry (with stochasticity-based weights, the effect should be larger for uniform 
weights such as the ones used by the authors) still tends to disperse clusters with real log(age) 
close to 7.5 into a distributi on significantly b roader than for clusters with log(age) of 7.0 or 
8.0. I suspect that Fig. 2 of iFall et al.l (120051 ) is biased by this effect and that a significant 
fraction of the clusters in the second age bin in reality belong to the third. 

Therefore, this paper shows that the existing analyses for M33, M51, and the Antennae 
likely introduce age biases that produce in turn systematic errors in the age-mass plane. As 
a consequence, the evolution of the CMF with age derived from studies that do not take 
into account the appropriate biases should be considered with caution. In particular, the 
region between 10 Ma and 100 Ma is especially subject to SIMF sampling effects. This is 
unfortunate for the sake of evaluating the validity of the Baltimore and Utrecht models, 
since tha t is the age r ange where both models differ more strongly in their predictions (see 
Fig. 1 in iLamerd l2008l ) . 



7.4. Possible solutions 

The reader of the previous subsections may get the impression that using unresolved 
photometry to obtain accurate ages and masses for stellar cluster systems is a hopeless 
endeavor. However, I believe that is not entirely the case if certain steps are taken. Some of 
the ideas below are derived directly from the results in this paper while others need further 



testing: 



1. Monte Carlo simulations (such as the ones in this paper) are necessary and 

should be performed for the specific filter sets used and ranges in mass and age studied. 
The results of the simulations are needed to estimate the mass and age ranges where 
results are valid and to generate Bayesian priors (see next point). 

2. An advanced Bayesian procedure should be applied to the likelihood mass- 
age map. Most analyses (including the one in this paper) simply obtain the mode of 
the likelihood in mass and age (i.e. the values that minimize x^) ^i^nd take those values 
as the observed ones. A further step would be to evaluate the likelihood of all possible 
values, determine a joint mass-age distribution using a reasonable Bayesian prior, derive 
a combined likelihood for all clusters (i.e. a mass-age distribution function), use the 
latter as a new Bayesian prior, and iterate as needed. In this way, it may be possible 
to eliminate some of the age biases discussed here, especially in those cases where 
multimodal solutions exist. 

3. Weights should be assigned in such a way that filters with lower stochasticity 
contribute more to the likelihood map than those with higher stochasticity. Note that 
in this paper I have used specific weights for each age, which is not possible with real 
data because we do not know the real ages to start with0. However, as it is clear from 
Table [3], the dependence on wavelength of 6 is rather similar for most ages (except for 
1 Ma clusters but, as we have seen, age fitting for those clusters is quite robust), so 
the use of age-averaged weights should not introduce large changes. 

4. Select the filters before observations are performed. One advantage of perform- 
ing the Monte Carlo simulations in advance is that it is possible to select the filters 
that maximize the information on the stellar content. In particular, 1 hypothesize 
(but do not prove) that the use of in termediate- and narr ow-band filters that measure 



quantities similar to the Lick indices (IWorthey et al.lll994j ) but not necessarily confined 



to that region of the spectrum (e.g. they could include Ha, other TiO bands, or the 
Ca triplet) could lower the values of 6, especially when age and extinction need to be 
disentangled. The use of such filters may be hard to implement because of the simulta- 
neous need for a large collecting area (compared to broad-band photometry) and high 
spatial resolution over large fields of view and wavelength coverage (which tends to 



""^^But it is not out of the question to use an iterative procedure fitting the observed ages such as the one 
in the previous point to address this issue. 
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favor space telescopes over ground-based techniques). Intermediate- and narrow-band 
filters are also helpful in combination with the idea in the next point. 

5. Incorporate the Monte Carlo simulations into the age-fitting code. This 
implies carefully selecting a relative large number of SIMF realizations (^100) for 
each age and mass range, computing the predicted photometry for all of them, and 
running a CHORIZOS-like code with the realizations as an added dimension. Then, 
after the likelihood is computed, the additional dimension can be collapsed in order to 
end up with a likelihood map (in e.g. mass, age, and extinction) that can be analyzed 
with Bayesian techniques. Such an analysis would effectively take into account SIMF 
sampling effects but is likely to be very costly in computational terms. The use of 
intermediate- and narrow-band filters would be a nice inclusion into such a code because 
the measurement of e.g. the Ca triplet equivalent width would allow an estimation of 
the real number of RSGs and AGB stars when the predicted number for their number 
at a given mass and age is so low that SIMF sampling effects are very large. 

6. If the clusters can be partially resolved into stars, use the number and colors 
of the bright point sources to con strain the giant /s uper giant population. 



See the analysis for NGC 4214 I-Es in lUbeda et al 



fl2007al jbr) "or an example. The 



objections to such an strategy are that it is labor intensive (it requires a one by one 
detailed effort for each cluster) and that with HST (or similar) resolution it can be 
extended effectively only up to distance s of a few Mpc. Nevert heless, it would be an 



interesting project to e.g. reprocess the IChandar et al.l (jl999bl ) results for M33 with 



such an strategy and adding information from longer wavelength data obtained with 
new observations of the same fields. 

7. Also, for partially resolved clusters and associations, attempt to do separate 
fits to the existing subclusters. As previously mentioned, multiple ages associated 
to different regions in extragalactic stellar clusters have been known to exist for some 
time now. Therefore, whenever high-spatial-resolution imaging is available, one should 
try to find age differences between subgroups for young objects. 



8. Conclusions 

1. SIMF sampling effects are the largest contributor to the observational uncertainties 
of the observed masses (and ages) of unresolved clusters studied with broad-band 
photometry. 



2. A Poisson mass (or equivalent quantity) needs to be calculated for each set of used 
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filters and possible values of the additional parameters (age, extinction. . . ) in order 
to determine the range of validity of the results for the individual cluster masses and 
the CMF. If multiple filters are used, the most straightforward way to calculate the 
Poisson mass is with Monte Carlo simulations. 

3. A rule of thumb is that if the real CMF slope is close to —2.0 and one desires to 
calculate the CMF slope with an accuracy of 0.1 without introducing any corrections, 
then only masses larger than 10 times the Poisson mass can be included. 

4. For most ages, the Poisson mass of a given filter increases with wavelength between U 
and K. 

5. The use of filter-dependent weights based on the variation of the Poisson mass with 
wavelength can reduce the biases and uncertainties introduced by SIMF sampling. 

6. When using broad-band optical+NIR photometry, it is possible to accurately measure 
masses and ages of young clusters (MS, LBV, and RSG phases) down to 10*^ Mq and 
possibly less. 

7. Post-RSG clusters with ages less than 100 Ma are especially affected by SIMF samphng 
effects. 

8. The behavior of SIMF sampling effects in old clusters is strongly dependent on whether 
extinction is known or not and on the possible existence of UV-bright stars (pAGB, 
HE, blue stragglers). 

9. There are several avenues of improvement that could reduce the uncertainties and 
biases induced by SIMF sampling. 

10. But one also has to recognize that in the low cluster-mass hmit there are cases where 
it is not possible to unambiguously measure individual cluster masses from unresolved 
photometry with accuracy because two clusters with different masses can have near- 
identical luminosities and colors. E.g. two clusters, one with a mass of 100 Mq and 
another one with a mass of 200 Mq, that only differ in the existence in the second one 
of a 100 Mq star (or of a 60 Mq star and a 40 Mq star) that has already exploded as 
a SN. In such limit, only an statistical interpretation of the CMF should be possible. 
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A. A Bayesian formalism for the compcirison between observed and real 

quantities 

Suppose that for each member of a population of astronomical objects wc want to 
measure a quantity X (luminosity, mass, distance. . . ) and study the properties of its real 
distribution fx{x). An object with a real value of X given by x will be observed to have 
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a value of Xo, with Xo ^ x in general. The observed quantity Xo may differ from the real 
one X because of observational or model uncertainties. If uncertainties can be described by 
a Gaussian with mean x and standard deviation a, then the distribution of Xo for a given x 
(i.e. the conditional distribution of the observed quantity) will be given by: 



fx^\x{Xo\x) 



2'Ka 



exp 



a 



(Al) 



the joint distribution of Xo and x (i.e. the global distribution in a 2-dimensional observed + 
real space) by: 



fx,,x{Xo,x) = fx,\xiXa\x)fxix) 



2Tia 



exp 



fxix), (A2) 



and the observed distribution by: 



fxSxo) = j fx,,x{xo,x) dx 



2na 



exp 



1 f Xo — X 
2 



fx{x)dx. (A3) 



The combination of [a] fxSxo) being different from fx{x) and [b] the possible existence 
of selection effects linked to e.g. s ensitivity limits is the source of a number of biases in 
astronomy, such a s the Lutz-Kelker Jhutz fc Kelkeil l 19731: ISmithll2003h and Malmquist biases 
(jMalmquistlll922l : iButkevich et al.ll2005l ). Biases affect the statistical properties of a sample 
of observed objects and they depend not only on the observed values and uncertainties 
(as it is sometimes as sume d) but also on selection effects and the properties of fx{x) (see 
Mafz Apellanig l2001al and ISmithI l2003l for the Lutz-Kelker case). If fx{x) is known, it is 
possible to compute a correction for individual objects using: 



X 



Xc 



J fxo,x{xo,x) xdx 
J fxo,x{xo,x) dx 



(A4) 



where I am allowing for the possibi lity that fxo,x{xo, x) is not normalized (see lMaiz Apellaniz et al 
2004bl : iMaiz Apellaniz et al.ll2008l for examples of applications of such corrections to trigono- 
metric parallaxes). Our best estimate for X would then be Xq + c. In practical terms, this is 
complicated by the fact that c itself depends on fx{x). One solution is to assume a param- 
eterized functional form for fx {x) with a reasonable initial gu ess for the parameters and t o 
iterate over the full procedure until convergence is reached (see iMaiz Apellanizll2001al . 120051 ) . 
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Besides knowing a correction for the value of an individual object, it is also interesting 
to know whether the observed distribution at a given value is larger or smaller than the real 
distribution i.e. whether uncertainties contribute to increasing or decreasing the number of 
objects at a certain point in x^. This can be evaluated through the fractional number change 
p between real and observed distributions, given by: 



fx{Xo) Jx{Xo) 

In this paper I apply the formalism above to the measurement of CMFs from photometric 
data. I will assume that the real cluster mass distribution is given by a power law: 



fM{m) = Arn^ (A6) 

between a lower and an upper mass limits miower and mapper, with A being a normalizing 
factor and 7 the power law exponent!^. 

In general, if fhiif^) is given by a power law, /a/oI^o) will be a different function. 
However, we can still define a local value for the slope 7' of fMoi^o) and a change in slope 
A7 between the real and observed distributions as: 

, dfM,{rno)/dmo 
A7 = 7 - 7 = — — — — 7. (A7) 



^^Note that some papers use /? for the exponent. Here I prefer 7 to maintain notation uniformity with the 
stellar case. 

Table 1. Mass uncertainty sources for the Antennae observations described in the text. 



Type Mass dependence cr (Mq) 

m = 10^ M„ m = 10^ M, 



Photometric m°-^ (approx.) 2.6 • 10^ 2.2 ■ 10^ 

Model complex < 5.0 ■ 10^ < 5.0 ■ 10^ 

SIMF sampling m"-^ 7.4 ■ 10^ 7.4 ■ 10^ 
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Table 2. Results of the Monte Carlo simulations of /mo|m('^o|'^) ioi 10 Ma clusters 
observed in V with known age and extinction. 
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s 


m 





a 




e 
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71 


7i,r 


72 


72,r 


(Mo) 






(M 


^) 


(Mo) 


(Mo) 














3 


■ 102 


3.3- 


10^ 


3.01 ■ 


102 


4.084 ■ 


102 


556.1 


5.302 ■ 


10-" 


3.694 


2.747 


17.603 


11.316 


1 


■ 10^ 


1.0 ■ 


10-^ 


1.00 ■ 


10^ 


7.494 ■ 


102 


561.6 


1.767- 


10+0 


2.022 


1.504 


5.238 


3.395 


3 


■103 


3.3 • 


10^ 


3.00 • 


10'^ 


1.300- 


103 


563.1 


5.302 - 


10+0 


1.186 


0.869 


1.785 


1.132 


1 


■10" 


1.0 ■ 


10" 


1.00 ■ 


10" 


2.397- 


10^ 


574.7 


1.767- 


10+1 


0.643 


0.476 


0.456 


0.339 


3 


■ 10^ 


1.0 ■ 


10" 


3.00 • 


10" 


4.133- 


10^ 


569.4 


5.302 - 


10+" 


0.366 


0.275 


0.361 


0.113 


1 


• 10^ 


1.0- 


10" 


1.00- 


10^ 


7.488 - 


10^ 


560.6 


1.767- 


10+2 


0.202 


0.150 


0.038 


0.034 


3 


• 10^ 


1.0- 


10" 


3.00- 


10^ 


1.308 - 


10^ 


570.1 


5.302 - 


10+2 


0.155 


0.087 


-0.031 


0.011 


1 


• 10^ 


1.0- 


10" 


1.00- 


10^ 


2.389 - 


10^ 


570.9 


1.767- 


10+3 


0.110 


0.048 


0.004 


0.003 



Table 3. Values of 6 for different log(ages) (in years) and filters for clusters with known 

age and extinction. 



(Mo) 



Filter 


6.0 


6.5 


7.0 


7.5 


8.0 


8.5 


9.0 


9.5 


10.0 


U 


235 


945 


160 


149 


53 


33 


23 


26 


42 


B 


200 


1314 


192 


305 


51 


32 


27 


52 


92 


V 


190 


1555 


566 


370 


89 


38 


51 


106 


163 


R 


186 


1722 


1034 


443 


151 


72 


109 


173 


242 


I 


178 


1949 


1585 


616 


273 


463 


476 


477 


581 


J 


157 


2 741 


2 240 


1170 


708 


3 598 


2 930 


2 702 


2 445 


H 


146 


3 081 


2476 


1640 


1281 


5184 


4 828 


4499 


4123 


K 


140 


3 333 


2 542 


1771 


1527 


5 717 


5 688 


5 512 


5 058 



Table 4. Values of mo,LLL for different log(ages) (in years) and filters for clusters with 

known age and extinction. 



"^o,lll(M0) 



Filter 


6.0 


6.5 


7.0 


7.5 


8.0 


8.5 


9.0 


9.5 


10.0 


U 


1093 


5 741 


1745 


754 


223 


96 


57 


81 


175 


B 


1018 


13477 


3265 


1164 


288 


103 


193 


315 


577 


V 


1000 


17509 


2 796 


1192 


464 


410 


551 


650 


1074 


R 


993 


20 302 


2 795 


1203 


998 


896 


961 


1058 


1524 


I 


978 


23 911 


2 801 


1954 


2 755 


2 643 


2 501 


2 825 


4189 


J 


930 


36 216 


4083 


3 894 


7328 


6 070 


8 796 


14406 


19417 


H 


904 


41098 


4570 


5129 


11822 


7386 


11734 


19 235 


26 623 


K 


889 


44401 


4809 


5 572 


13 891 


7806 


13369 


21949 


30 617 
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Table 5. Ages and masses calculated with CHORIZOS for 10^ Mq clusters with unknown 
age and known extinction. The last column gives the expected mass uncertainty for known 

age and extinction. 



Real log(age/l a) m (Mq) a (Mq) 

log(age/l a) 7 filt. 7 filt. 3 filt. 7 filt. 7 filt. 3 filt. V 

unweig. unweig. 



6.0 


6, 


.25 


+0. 

-0, 


,08 
,17 


6, 


.25 


+0, 
-0 


,08 
.17 


6. 


.25 


+0, 
-0 


,08 
.19 


8765 


+1002 
-3204 


8753 


+ 1013 
-3164 


8657 


+ 1119 
-3085 


1378 


6.5 


6, 


.45 


+0. 

-0, 


06 

,05 


6, 


.45 


+0, 
-0 


06 
.05 


6, 


.45 


+0, 
-0 


,06 
.05 


9846 


+3537 
-2418 


9855 


+3527 
-2420 


9773 


+3631 
-2548 


3943 


7.0 


7, 


.05 


+0. 

-0, 


,04 
,15 


7, 


.05 


+0, 
-0, 


,07 
.04 


6. 


.93 


+0, 
-0 


12 
.34 


9558 


+2003 
-1789 


11179 


+3039 
-2573 


7879 


+3311 

-5575 


2379 


7.5 


7, 


.37 


+0. 

-0, 


,32 
,25 


7, 


.29 


+0, 
-0 


,39 
.18 


7. 


.47 


+0, 
-0 


39 
.26 


6514 


+7009 

-2549 


5903 


+7393 
-1983 


9427 


+9458 
-4960 


1922 


8.0 


7, 


.97 


+0. 

-0, 


14 
,22 


7, 


.90 


+ 0, 

-1 


19 
.21 


7. 


.99 


+0, 
-0 


,08 
.14 


9404 


+3163 
-3166 


8558 


+3712 
-7910 


9782 


+ 2285 
-2362 


945 


8.5 


8, 


.48 


+0. 

-0, 


,07 
,12 


8, 


.48 


+ 0, 
-0 


20 
.48 


8. 


.49 


+0, 
-0 


,04 
.03 


9604 


+ 1676 
-1671 


9543 


+4267 
-5021 


9951 


+ 1103 
-930 


620 


9.0 


8, 


.99 


+0. 

-0, 


,09 
,08 


9, 


.00 


+ 0, 
-2 


25 
.18 


8. 


.99 


+ 0, 
-0 


()5 
.05 


9789 


+2537 
-1748 


10194 


+6189 
-10027 


9852 


+ 1962 
-1413 


715 


9.5 


9, 


.48 


+0. 

-0, 


16 
,10 


9, 


.50 


+ 0, 
-0 


25 
.50 


9. 


.49 


+0, 
-0 


11 
.08 


9720 


+4494 
-2784 


9756 


+6754 
-6631 


9844 


+3594 
-2368 


1032 


10.0 


9, 


.95 


+0. 
-0, 


,05 
.17 


9. 


.94 


+0, 
-0, 


,05 
.44 


9. 


.98 


+0, 
-0, 


,02 
.12 


9003 


+2172 
-3428 


8632 


+2522 
-5323 


9413 


+ 1792 
-2869 


1277 
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Table 6. Ages, color excesses, and masses calculated with CHORIZOS for 10^ Mq clusters 
with unknown age and extinction. The last column gives the expected mass uncertainty for 

known age and extinction. 



Real log(age/l a) E(4405 - 5495) m (Mq) a (Mg) 

log(age/l a) 7 filt. 3 filt. 7 filt. 3 filt. 7 filt. 3 filt. V 



6.0 


6, 


.281°: 


14 

,28 


6, 


,35 


+0, 

-0, 


,40 
,16 


1, 


,00 


+0, 

-0, 


,00 

,04 


0, 


,99 


+0, 

-0, 


,01 
,54 


8542 


+ 1016 

-3466 


6806 


+2894 
-2142 


1378 


6.5 


6, 


,44 1«: 


,06 


6, 


,47 


+0, 
-0 


,37 
,06 


1, 


,00 


+0, 

-0 


04 
,01 


0, 


,99 


+0, 
-0 


,02 
,53 


10142 


+4720 
-2711 


11031 


+8596 

-3448 


3943 


7.0 


7, 


,051°: 


,01 
,11 


6, 


,90 


+0, 
-0 


,08 
,43 


0, 


,99 


+0, 
-0, 


12 
,15 


0, 


,98 


+0, 
-0 


,32 
,15 


10101 


+2393 
-2494 


6380 


+2653 
-3374 


2379 


7.5 


7, 


,471°: 


,25 
,34 


7, 


,47 


+0, 
-0 


,48 
,27 


1, 


,01 


+0, 
-0, 


14 

,11 


1, 


,00 


+0, 
-0 


,10 
,10 


8774 


+6047 
-4335 


9319 


+ 11052 
-4953 


1922 


8.0 


8, 


■001°: 


,00 
,21 


7, 


,99 


+0, 
-0 


,05 
,09 


1, 


,00 


+0, 

-0, 


,08 
,07 


1, 


,00 


+0, 
-0 


,06 
,06 


9461 


+ 1203 
-1939 


9569 


+854 
-910 


945 


8.5 


8, 


,53 IJ: 


13 

,03 


8, 


,85 


+0, 
-0 


,93 
,33 


0, 


,96 


+0, 

-0 


09 
,66 


0, 


,76 


+0, 
-0 


,22 
,50 


10807 


+83593 
-1203 


18376 


+ 103069 
-7893 


620 


9.0 


9, 


,34 «: 


,39 
,90 


9, 


,16 


+0, 
-0 


,33 
,72 


0, 


,80 


+ 0, 
-0 


59 
,23 


0, 


,94 


+0, 
-0 


,41 
,25 


20898 


+22283 
-17097 


15029 


+ 12787 
-11162 


715 


9.5 


9, 


.771°: 


,19 
,66 


8, 


,94 


+ 0, 
-0 


,57 
,44 


0, 


,90 


+0, 

-0, 


30 
,11 


1, 


,34 


+ 0, 
-0 


,31 
,35 


16520 


+6497 
-10980 


3170 


+6936 
-1730 


1032 


10.0 


10, 


.001°: 


,00 
,20 


9, 


,98 


+0, 
-0, 


,02 
,04 


0, 


.97 


+0, 
-0, 


14 
,07 


1, 


,00 


+0, 
-0, 


,05 
,05 


9305 


+ 1497 
-2361 


9251 


+ 1411 
-1241 


1277 
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Fig. 1. — Results of the Monte Carlo simulations of /mo|m(^o|^) and Gaussian and Gamma 
fits for 10 Ma clusters observed in V with four real masses: 300 Mq (upper left), 1000 Mq 
(upper right), 10 000 M© (lower left), and 100 000 M© (lower right). The vertical dotted 
lines show the values of the real mass m and, when visible, of 9 and moLLL- 
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Fig. 2. — Observed [/moI^^o), continuous lines] and real [/Af(m), dotted lines] mass distri- 
butions (top) and slopes (7', bottom) using the Gamma approximation with 10 Ma clusters 
observed with V {6 = 566 Mq). 
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Fig. 3. — Isochrones for seven of the nine ages used in this paper. For each isochrone, 
symbols of different colors are used to indicate the median mass for the luminosity in the 
eight filters U, B, V, R, I, J, H, and K assuming a Kroupa IMF between 0.1 Mq and 120 
Mq. The two dashed lines join the points in each isochrone beyond which 0.20% and 0.10%, 
respectively, of the remaining stars in a well-sampled Kroupa IMF are located. The two 
remaining isochrones are not included for the sake of clarity. 
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Fig. 4. — m*,av/^eff and mo,LLL as a function of 6 for the 72 Monte Carlo simulations (9 ages 
times 8 filters). The straight line marks the 1 to 1 locus. 
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Fig. 5. — Observed [/moI^^o), continuous lines] and real [fM{fn), dotted lines] mass distribu- 
tions using the Gamma approximation with [a] 1 Ga clusters observed with B {6 = 23 Mq, 
top) and [b] 316 Ma clusters observed with K [9 = 5717 M©, bottom). 
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Fig. 6. — Distribution of observed ages for the nine ages and three execution types described 
in the text for clusters of 10^ Mq and known extinction. At the top of each panel the dashed 
line indicates the real age while the symbols and error bars provide the median and inferior 
and superior uncertainties (1-sigma equivalents) of each distribution. 
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Fig. 7. — Same as Fig. [6] for clusters of 10^ Mq and unknown extinction but only for the 
first and last execution types. 
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Fig. 8. — Same as Fig. [6] for clusters of 10^ Mq and unknown extinction but only for the 
first and last execution types. 



